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ABSTRACT 

We examine radial oscillations of superfluid neutron stars at finite internal 
temperatures. For this purpose we generalize the description of relativistic 
superfluid hydrodynamics to the case of superfluid mixtures. We show that 
in a neutron star at hydrostatic and beta-equilibrium the red-shifted temper- 
ature gradient is smoothed out by neutron superfluidity (but not by proton 
superfluidity). We calculate radial oscillation modes of neutron stars assuming 
"frozen" nuclear composition in the pulsating matter. The resulting pulsation 
frequencies show a strong temperature dependence in the temperature range 
(0.1 — l)Tcn, where Ten is the critical temperature of neutron superfluidity. 
Combining our results with thermal evolution, we obtain a significant evolu- 
tion of the pulsation spectrum, associated with highly efficient Cooper pairing 
neutrino emission, for 20 years after superfluidity onset. 
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1 INTRODUCTION 

It is commonly accepted that a neutron star becomes superfluid (superconducting) at a 
certain stage of its thermal evolution (see, e.g., Lombardo and Schulze 2001). It is believed, 
in particular, that protons pair in the spin singlet (^5*0) state, while neutrons pair in the 
spin triplet ('^^2) state in the neutron star core. A large number of different models of 
nucleon pairing have been proposed in literature (references to original papers can be found 
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in Yakovlev et al. 1999 and in Lombardo and Schulze 2001). These models predict very 
different density profiles of neutron (n) and proton (p) critical temperatures, Tcn(p) and 
TcpIp), respectively. 

In spite of the many theoretical uncertainties in the theory, it is clear that superfluidity 
strongly affects the neutron star evolution, for example, its cooling (see, e.g., Yakovlev and 
Pethick 2004, Page et al. 2004), neutron star pulsations (sec e.g., Mendell 1991a,b, Lindblom 
and Mendell 1994, Lee 1995, Andersson and Comer 2001a, Andersson et al. 2002, Prix et 
al. 2004), and is probably related to pulsar glitches (see Alpar et al. 1984, Andersson et al. 
2003, Mastrano and Melatos 2005, Peralta et al. 2005). 

In this paper we discuss the effect of superfluidity on neutron star dynamics. The hydro- 
dynamics of a superfluid liquid, composed of identical particles, was formulated by Khalat- 
nikov (1952) within Tisza's (1938) two-fluid model, which was elaborated by Landau (1941, 
1947). This "orthodox" two-fluid model is based on the assumption of two independent ve- 
locity fields: the "normal" velocity of thermal excitations Vq and the "superfluid" velocity 
Vg, each carrying some part of the mass of liquid, so that the mass current density j can be 
written as 

3^{p- P.)V^ + p,V,, (1) 

where ps is known as the superfluid density. The superfluid component moves without fric- 
tion and does not interact with the normal fluid. The hydrodynamic equations in this case 
include the equation of motion for the superfluid component, in addition to the energy and 
momentum conservation laws and the continuity equations for mass density and entropy 
(see, e.g., Putterman 1974, Landau and Lifshitz 1987, Khalatnikov 1989). 

Obviously, the hydrodynamics described above cannot be applied directly to superfluid 
neutron stars. The stellar core consists of, at least, three kinds of particles (neutrons, protons, 
and electrons), and neutrons and protons may be superfluid. The superfluid hydrodynamics 
was extended to superfluid mixtures by Arkhipov and Khalatnikov (1957) and Khalatnikov 
(1973) and later, more accurately, by Andreev and Bashkin (1975). 

The main element of hydrodynamics and kinetics of superfluid mixtures is the entrain- 
ment matrix pj^, which naturally appears in the theory as a generalization of the superfluid 
density to the case of superfluid mixtures. If the only baryons in the core are neutrons 
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and protons, the matrix p^k can be found from the relations (Andreev and Bashkin 1975): 

in = (Pn - Pnn " Pnp) V'q + Pnn V'ns + Pnp V"ps , (2) 

ip = (Pp - PPP - Ppn) + Ppp ^ps + Ppn t^ns • (3) 

Here pi — miUi, rrii is the mass of a free particle and rij is the number density of particle 
species i with i = n or p; and Vjg are the mass current density and the superfluid velocity of 
particle species i, respectively. Since "normal" protons and neutrons will be locked together 
by friction, we assume that their velocities Vq are identical. In other words we assume 
that the characteristic time Tnp of neutron-proton collisions is negligible in comparison with 
the typical hydrodynamic time (e.g., the inverse frequency uj~^ of stellar pulsations). For 
example, for non-superfluid matter Tnp ~ (10~^^ — 10~^^) Tg^ s (see, e.g., Yakovlev and 
Shalybkov 1991) is much smaller than u;~^ ~ 10~^ s (see Section 6), where Tg — T/{10^K) 
is the temperature in units of 10^ K. 

It follows from the phenomenological analysis of Andreev and Bashkin (1975) that the 
matrix pik is symmetric: pnp = Ppn- Moreover, at zero temperature the equalities 

Pnn + Pnp = Pn, Ppp + Ppn = Pp (4) 

must hold (see, e.g., Borumand et al. 1996), in order for the system to be invariant under 
Galilean transformations. 

The cntrainment matrix pik for a non-relativistic neutron-proton mixture was calculated 
by Borumand et al. (1996) at zero temperature and by Gusakov and Haensel (2005) for any 
temperature. At T = entrainment coefficients analogous to the matrix pi^ have also been 
calculated by Comer and Joynt (2003). Even though neutrons (and certainly protons) can be 
considered non-relativistic with good accuracy up to the densities p < 10^^ g cm~^, a fully 
relativistic calculation of Comer and Joynt (2003) is more self-consistent. Nevertheless, we 
will use the results obtained by Gusakov and Haensel (2005), because we deal with dynamic 
effects associated with finite temperatures in the neutron star core. 

The hydrodynamics of superfluid mixtures presented by Andreev and Bashkin (1975) 
cannot be applied directly to neutron stars, since it is an essentially non-relativistic theory. 
We need to generahze the description to take into account the effects of General Relativity 
which are important to neutron stars. Landau's two-fluid model, initially applied to liq- 
uid helium 11, was extended to General Relativity by Carter (1976, 1979, 1985) using a 
convective variational principle and by Khalatnikov and Lebedev (1982) and Lebedev and 
Khalatnikov (1982) on the basis of a potential variational principle. The equivalence of these 
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two approaches in the non-dissipative hmit has been demonstrated by Carter and Khalat- 
nikov (1992a,b). The former approach was extended by Carter and collaborators to analyze 
superfluid mixtures, in particular, neutron star matter (see, e.g., Langlois et al. 1998, Carter 
and Langlois 1998). The hydrodynamic equations derived from the convective variational 
principle impose restrictions on the "canonical" coordinates and momenta. They include 
various phenomenological coefficients, which need to be related to parameters that are actu- 
ally calculated from microscopic theory, for example, the superfluid densities. For this reason 
we will not use Carter's elegant framework here, even though all available calculations of 
superfluid oscillations in General Relativity have so far been made within this approach (see 
Comer et al. 1999, Andersson and Comer 2001b, Andersson et al. 2002, Yoshida and Lee 
2003). Instead, we will employ a version of superfluid hydrodynamics derived by Son (2001) 
from microscopic theory (see also Pujol and Davesne 2003, Zhang 2003). Slightly modified, 
this approach has the advantage of offering an easy interpretation of the various physical 
quantities entering the hydrodynamic equations. Although one can show that our equations 
are formally equivalent to those of Carter, it is clear that further work is needed to connect 
his formulation with the microphysics. 

The aim of the present study is to analyze the effect of finite temperatures on pulsations 
of superfluid neutron stars. Pulsations may be excited during the star's formation or during 
its evolution under the action of external perturbations (e.g., accretion, gravitational per- 
turbations) or internal instabilities (associated with unstable pulsation modes). A possible 
signature of these pulsations would be the modulation of the electromagnetic radiation from 
the neutron star surface or the detection (in the future) of gravitational radiation generated 
by nonaxisymmetric fluid motion. It will be shown that the effect of finite temperatures 
may essentially influence the pulsation spectrum in the temperature range T ~ (0.1 — l)Tcn, 
because in this range the entrainment matrix pik changes considerably and cannot be treated 
as a constant. This, in turn, affects the hydrodynamic equations for superfluid mixtures and, 
hence, the oscillations of the star. To simplify the problem, we restrict ourselves to the case 
of radial pulsations and examine a simple one-fluid model of the non-elastic neutron star 
crust consisting of normal matter. The core will be assumed to consist of neutrons, protons 
and electrons (npe- matter) , with both types of nucleons being superfluid. 

We would like to note that all previous calculations of global pulsations of superfluid 
neutron stars were made in a zero-temperature approximation. We believe that this is too 
idealized for two reasons. First, even an initially cold star can be heated by pulsations 
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because of the transformation of the pulsation energy into heat (for example, due to viscous 
dissipation, see Gusakov et al. 2005). Second, the critical temperatures of nucleons depend 
on the density. This is a bell-shaped curve, which shows that the critical temperature first 
rises with the density and then decreases after reaching a maximum. Thus, for any given 
temperature T there is usually a region in the star with T ~ Ten- This is an important point 
that is worth emphasizing. 

This paper is organized as follows. In Section 2 we extend Son's equations to super- 
fluid mixtures and rewrite them using more appropriate variables. In Section 3 we consider 
equilibrium configurations of neutron stars. In Section 4 we discuss equations for radial pul- 
sations taking into account a finite temperature in the core. In Section 5 we analyze short 
wavelength solutions to these equations, i.e., sound waves in the superfluid neutron star. In 
Section 6 we examine the numerical solutions to the pulsation equations and the eigenfre- 
quency spectrum as a function of temperature. In addition, we study the evolution of the 
oscillation spectrum during the star cooling. 

2 RELATIVISTIC EQUATIONS FOR NON-DISSIPATIVE 
HYDRODYNAMICS OF SUPERFLUID MIXTURES 

In this section, the relativistic equations suggested by Son (2001) for a one-component 
superfluid liquid at finite temperature will be extended to multicomponent mixtures, and 
rewritten in a form which is better suited for our application. For simplicity, let us consider 
a mixture of three kinds of particles, assuming that two kinds are superfluid and one kind 
is normal. In a neutron star, for example, neutrons and/or protons may be superfluid, while 
electrons (with species index e) are normal. 

It is well known that, in superfluid matter, several independent motions with different 
velocities may coexist without dissipation (see, e.g., Khalatnikov 1989). When a mixture is 
composed of two superfluids and one normal fluid (in principle, there may be many normal 
species), the system is fully deflned by three 4- velocities u^, '^{■a)^ ^""^ '"^fp)- '^^^ latter two 
arise from additional degrees of freedom associated with superfluidity. The velocity refers 
to electrons as well as "normal" neutrons and protons (Bogoliubov excitations of neutrons 
and protons). 

If there are several independent motions, the question arises how to define the comoving 
frame in order to determine the basic thermodynamic quantities: the energy density e and 
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the particle number densities rii {I = n, p,e). Without any loss of generality, we can assume 
that the reference frame in which the velocity u'^ equals to u'^ = (1, 0, 0, 0) is comoving. 
This assumption imposes certain restrictions on the particle 4-current and the energy- 
momentum tensor T'^'^ 

= UimUuT'"' = e. (5) 

The full set of hydrodynamic equations for superfiuid mixtures which satisfy these conditions 
is 

de = TdS + drii + /Xe dng + ^ d (w(l)W(fc)a) , (6) 
J'?);^ = 0' j[i)=niu'' + YikW^^^^, j^^^ = n,u^', (7) 
^ T'^'^ = {P + e)u'^u-' + Pg'^'' + Yik[w^^wl,)+^i,w^,^u' + ^i,w'(,)u''), (8) 
u,w^,^ = 0. (9) 

Here and below, the subscripts i and k refer to nucleons: i,k = n, p. Unless otherwise stated, 
a summation is assumed over repeated spacetime indices (Greek letters) /x, z/, a and nucleon 
species indices (Latin letters) i, k. Eq. (jH)) represents the second law of thermodynamics 
for superfiuid mixtures, while Eqs. ((Zj) and (jH)) describe particle and energy-momentum 
conservation laws, respectively. Finally, Eq. © is the additional equation for a superfiuid 
component; it is a necessary condition for Eq. © to hold. 

In Eqs. (0)-© 9^'^ is the metric tensor; S is the entropy per unit volume; /i/ is the 
relativistic chemical potential of particle species / = n, p, e; P is the pressure which is 
defined in the same way as for ordinary (non-superfiuid) matter: 

P = —e + HiUi + /ie"-c + ^"S*. (10) 

Finally, Yn^ = Y^i is a 2 x 2 symmetric matrix, whose elements are the functions of temper- 
ature T and the number densities of neutrons and protons. Using Eqs. ^ and we can 
write the Gibbs-Duhem relation for a superfiuid mixture: 

dP = SdT + ni dfii + ne d/Xe - d {w"i)W(k)a) ■ (H) 

The requirement of constant total entropy of the mixture imposes an additional con- 
straint on the 4- velocities w'^-y Namely, we obtain the correct hydrodynamic equations for 
a perfect superfiuid mixture if the 4-velocities w'^--^ have the form 

w^^ = ^ - - (12) 
where (pi is an arbitrary scalar function, A^^ is the 4-potential of the electromagnetic field. 
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and Qi is the electric charge of nucleon species i. It is easy to demonstrate that with the 
quantity w^.^ given by Eq. (fT^ . the set of equations leads to entropy conservation: 

{Su^\, = 0. (13) 

Obviously, the entropy is carried with the same velocity as the normal fluid, i.e., the 
entropy of the superfluid fraction in the mixture is zero. 

Let us now specify the physical meaning of the quantities 0j, u^, and Yik- For this aim 
we will examine how they are related in the non-relativistic limit to the superfluid velocity 
Vis, the normal velocity Vq, the wave function phase of the Cooper-pair condensate $i, and 
the entrainment matrix pik (these quantities appear in the non-relativistic hydrodynamics of 
superfluid mixtures discussed in detail by Andreev and Bashkin 1975). One can demonstrate 
that the following relations hold: 

1 h 

Vi, = — {V(t)i - QiA) , V</., = - V$i, (14) 

Vq = u, Yiu = Yki = -^. (15) 

Here and below, the speed of light is assumed to be c = 1. For convenience, a brief glossary 
of symbols is presented in Table 1. 

Let us discuss in more detail the properties of the matrix 1^^. In the absence of super- 
fluidity, when the temperature T is higher than the critical temperatures of neutrons T^n 
and protons Tcp, we have Yik = 0. Then the expressions for the 4-currents (|7j) and for the 
energy-momentum tensor (jHJ take the standard form and describe a normal perfect fluid 
(see, e.g.. Landau and Lifshitz 1987). If, for example, the inequality Ten < T < Tcp holds, 
i.e., if we have only superfluid protons, the only non-vanishing matrix element is Y^p. In 
contrast, at T = a// neutrons and protons form Cooper pairs. In other words, there are no 
nucleons moving with the normal fluid component at velocity u'^. A 4-current j^y therefore, 
is independent of u^, and we have the condition (see Eqs. 171 and IT^ : 

p,Yik{T = 0) = m. (16) 

Unfortunately, to our best knowledge, results for the matrix Yik at finite temperatures have 
not yet been presented in the literature. Nevertheless, Gusakov and Haensel (2005) calculated 
the entrainment matrix pik{T). As we have already mentioned, the matrices Yik and pik are 
interrelated by Eq. (fT3j) in the non-relativistic limit. Thus, we will use an approximate 
expression for the matrix Yik, which satisfies Eq. (fT3j) in the non-relativistic limit and at the 
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Table 1. A brief glossary of symbols describing superfiuid hydrodynamics in both the non-relativistic and the relativistic 
domains. Subscripts i and k refer to nucleons: i, k=n, p. 





critical temperature 


P 


density 


Ps 


superfiuid density 


Vq 


velocity of thermal excitations (Bogoliubov quasiparticles) 


Vs 


superfiuid velocity 


3 


mass current density 


T 


critical temperature of particles i 


Pi 


density of particles i 


Pik 


entrainment matrix 


*i 


wave function phase of the Cooper-pair condensate of particles i 




superfiuid velocity of particles i 


Ji 


mass current density of particles i 


rii 


number density of particles i 


Y,k 


relativistic entrainment matrix, in the non-relativistic limit = pik/fni'mk 




4- velocity of electrons and neutron and proton thermal excitations 


4>i 


scalar potential related to <E>i by V0i = h V3'i/2 


P- 

w,.., 


4- velocity which reduces to = r?ii(Vis — Vq) in the non-relativistic limit 


I** 


4-current of particles i 




4-potential of the electromagnetic field 


H 


electric charge of particles i 



same time meets the condition |T6|) : 

-i^ _y _ Pnp y _ Pnn "I" Pnp ~ f^n^^p^np yr _ Ppp + Ppn ~ f^pf^ii^pn /'lyN 

-'np -'^pn 1 ^ rm j-'pp 

The condition ()16|) can be derived from these formulas, if we take into account that Eqs. Q 
must hold at T = 0. 



3 EQUILIBRIUM CONFIGURATIONS OF SUPERFLUID NEUTRON 
STARS 

Let us now use the above formulas to describe neutron stars. For simplicity, consider a non- 
rotating star. We will often refer to the results of the pioneering work of Chandrasekhar 
(1964) devoted to radial pulsations of non-superfluid stars in General Relativity. The metric 
for a spherically symmetric star, which experiences radial pulsations, can be written as (see, 
e.g., Chandrasekhar 1964) 

ds^ = -e^'dt^ + r2d^]2 + dr^, (18) 
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where r and t are the radial and time coordinates, respectively; dfi is a solid angle element in 
a spherical frame with the origin at the stellar center. The metric functions v and A depend 
only on r and t. The quantities referring to a star in hydrostatic equilibrium will be marked 
with the subscript "0"; in particular, the metric coefficients of an unperturbed star will be 
denoted as VQ{r) and Ao(r). 

In the equilibrium neutron star the measurable physical quantities (e.g., the number 
densities) must be time- independent. Thus, the continuity equation for electrons ((Tj) and 
the expression for the 4-velocity of the normal component 

yield (for a spherically symmetric star!) 

^0 = e-^o/^ u^ = u' = u^ = 0. (20) 
Next, the continuity equations for neutrons and protons ((Tj) give 

= = = 0- (21) 
Finally, in view of Eq. (j^ . one obtains from Eq. (0) 

wf^ = 0. (22) 

It is clear from Eqs. ((201)^(1221) that the energy-momentum tensor (jHl) of an equilibrium 
superfluid star is the same as that of a non-superfluid one. Therefore, the formulas that 
describe hydrostatic equilibrium of non-superfluid stars can be applied to our case as well. 
In particular, the following formula is valid (see, e.g., equation 21 of Chandrasekhar 1964) 

New information can be obtained from Eq. ((221) • When written for neutrons, it gives, together 
with Eq. dni), 

^ = -/inoe-/^. (24) 



dt 

On the other hand, from Eqs. (fT^ . ()20|) and ()21|) one flnds 



^nO 



0. (25) 



dr 

It follows from Eqs. ^ and ^ that 

A(/.„oe-/^)=0. (26) 

It should be emphasized that the application of conditions (21) and (22) to protons will 
not yield a constraint similar to Eq. (pUj) for /ipo, because Eq. (fT^ for the protons depends. 
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additionally, on the 4-potential of the electromagnetic field. We are not interested here in 
the relation between A^^ and /ipo that can be derived from Eqs. (j2II) and (j22I)- 

Assuming that a star at hydrostatic equilibrium meets, in addition, the quasineutrality 
condition, rZeo = ''^pO; one gets from Eqs. (fTUj) and (fTTj) 

Pq + So = finQflhO + + ^O'S'o, (27) 

-3— = ^bo \- rieo \- So — — , (28) 

dr dr dr dr 

where nbo = ""-no + '^po is the baryon number density; (5/io = /ipo + /ieo — /^no- Substituting 

the expression for duo/dr from Eq. into Eq. and using Eq. (|77j). one obtains 

dPo d/ino 1 / r I c ^ '^'^0 ^on^ 

^ = nto ^ - 2 i^f^orieo + Tq^o) (29) 

A comparison of Eqs. and leads to the equahty 

ueo ^ (5/^oe''«/^) + 5o ^ (Toe'^o/^) = 0. (30) 

Note that, to derive this formula we have considered a star in hydrostatic equilibrium (but 
not necessarily in thermal, diffusive, or beta-equilibrium). If we assume, in addition, that in 
some region of the star (i) the thermal equilibrium condition is fulfilled 

i (Toe-/^) = 0, (31) 
and (a) neutrons are superfiuid, then Eq. ()30|) tells us that this region must be in diffusive 
equilibrium (the opposite statement is also correct: diffusive equilibrium means thermal 
equilibrium for the problem in question). Indeed, in this case we have from Eqs. pHjl and 

^ (/ino e'"/') = 0, ^ [(/ipo + /ieo) e^'^'] = 0. (32) 
These conditions describe the diffusive equilibrium of npe-matter and are quite standard 
(see, e.g.. Landau and Lifshitz 1980). The second condition of Eq. (jHSI) is nothing but a 
sum of the diffusive equilibrium conditions written for protons and electrons. Each of them 
includes a self-consistent electrostatic potential to ensure quasineutrality (the most recent 
discussion of diffusive equilibrium as applied to npe-matter of neutron stars is given by 
Reisenegger et al. 2006). We are not interested here in determining this potential: it cancels 
out after the summation. 

In this paper we assume that an unperturbed star is at hydrostatic and beta-equilibrium 
(i.e. (5/io = 0). In this special case one immediately obtains from Eq. ()30|) the thermal 
equilibrium condition ()31|) . Thus, we arrive at the conclusion that a (red-shifted) temperature 
gradient cannot exist in any region of a hydrostatically and beta-equilibrated neutron star 
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which contains superfluid neutrons. This situation is identical to that for pure helium II 
(see, e.g., Khalatnikov 1989). Note that, proton superfluidity imposes no such restrictions 
on the temperature gradient. 



4 RADIAL PULSATIONS OF SUPERFLUID NEUTRON STARS 

In this section we consider a star with small radial perturbations. Accordingly, in all the 
equations we shall neglect the quantities which are second order and higher in the pulsation 
amplitude and retain the linear terms. In addition, we will use the hypothesis of a frozen 
nuclear composition, neglecting the effect of beta-processes on the chemical composition of 
the core during the pulsations. This assumption is justified if the radial pulsation frequencies 
are uj ^ 1/r, where r is the characteristic time of beta-equilibration. Recall that for non- 
superfluid matter and under the condition |/Up + /le — /Un| ^ T it can be estimated that 
r ~ Tg"^ months if beta- relaxation proceeds via the modified Urea process (see, e.g., Yakovlev 
et al. 2001); for superfluid matter beta-relaxation rates were calculated by Haensel et al. 
(2000, 2001), and by Villain and Haensel (2005). The final assumption we make is the 
validity of the quasineutrality condition in a pulsating star, 

rie = Up, (33) 

which should hold since u is much smaller than the plasma frequency of electrons, Upe- In 
the following, the quantities containing no "0" subscript refer to a perturbed star. If A is 
a physical quantity in a perturbed star and Aq the same quantity in the unperturbed star, 
then we denote A — Aq = 6 A. 

The quasineutrality condition leads to equal 4-currents of electrons and protons: 

By substituting the expressions for the currents from Eq. 0, one gets 

YpkW^,^ = 0. (35) 

We will also need the continuity equation for baryons, which can be found by summing the 
continuity equations ((Zj) for protons and neutrons. With Eq. (jSSI), we obtain 

(nbU^ + Fn^t^f,)) =0. (36) 
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4.1 Basic equations 

Using the metric (fTHj) . one can write the linearized 4- velocity as 

v!' = e-'''\ M^ = e-^«/2^, ^2 = ^=^ = 0, (37) 

where v = dr/dt is the velocity of the normal component of the mixture in the radial 
direction. (Note a misprint in formula 25 of Chandrasekhar 1964: the expressions for m° and 
Uq must have u instead of uq.) Using Eq. (jSTj), one can find directly from Eq. (jU)) 

= 0. (38) 

In addition, because particles move only in the radial direction, we have 

wf^ = wf,^ = 0. (39) 

Therefore, the only non-zero components of the energy-momentum tensor are 

To° = -e, T,' = Ti = Ti = P, (40) 

= -{Po + 6o)v + AT,\ (41) 
rpo ^ e^"-'''' {Po + eo) V ~ e^"-"'' ATq\ (42) 

These formulas differ from those for the normal liquid only by the term ATq^, which is (see 
Eq.ED 

ATq = fikoYik uowli) = -finoYni e""^"^. (43) 

When writing the last equality, we have used Eq. and the expression uq = —e"/'^ (Eq. 1^1 
yields w^-^ ~ v, so that u can be substituted for z/q in Eq. l43|) . Note an important consequence 
of Eq. pnj) : if neutrons in a star are normal (Fni = 0), its pulsations will be indiscernible 
from those of a common non-superfluid star, no matter whether the protons are superfiuid 
or not. 

Let us analyze Eq. ()38p for neutrons. With Eq. ()12p it can be rewritten as 

-e-'^^-/Xne-'^/^ = 0. (44) 

By substituting z/ = z/q + (5z/(r, t), 0n = 0no + S(pn{f^, t), fJ'n = /^no + Sf^ni^, t) into Eq. (jH)) and 
using Eq. (jUj), we get 

-(5/in + \^^n^5u\ e'^"/^ (45) 



dt V 2 

On the other hand, in the linear approximation and in view of Eq. (j25|) . we have 

^(n) = e -Q^-^^nU =e - /inoe V. (46) 



By combining Eqs. (jl^ and 
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, we find 



(47) 



dt ^ ^'^^ ' ^ dr 

Let us introduce new variables Zi and ^ according to (there is no summation over i here!): 

dzi 



w 



(0 



dt 



dt 



(48) 
(49) 



The time integration of Eq. gives the relation between the variables and z-p 

limyx>kZk = Q- (50) 
Assuming now that all perturbations vary with time as exp(zcijt), we rewrite Eq. (jTrf) in the 



form: 



A^no e 



Xo-vii/2 



d_ 

dr 



uo/2 



1 



V 



(51) 



Thus, we have derived one of the equations that describe pulsations of a relativistic superfluid 
star. There is no analogue of this equation for non-superfluid stars. In order to determine 
the unknown eigenf unctions Zi and ^ and the frequency spectrum, it is necessary to find an 
additional pulsation equation. In principle, this can be done by writing Einstein's equations 
with the energy-momentum tensor given by Eqs. (001)^(021) • However, the situation can be 
considerably simplified because this energy-momentum tensor does not essentially differ from 
that used by Chandrasekhar (1964) in the analysis of pulsations of non-superfiuid stars (see 
his Eqs. 27 and 28). By adjusting his derivation to our case, we find the following expressions 
for the quantities 5A, fe, and d5v/dr: 



5X 

6e 
d6iy 



1 d ,^ 



° Po + eo dr 
r"^ dr \ / ' 



dr r 



dr Po + eo 

Here, the quantity is defined by 



^ (Ao + i^o) ■ 



-^0 



dTl 
dt 



(52) 
(53) 
(54) 

(55) 



and is found to be (see Eqs. EJ 1121 UHl and | 

= — (Pq + ^o) ~ t^nOf^iO^ni Zi- (56) 

Eqs. are generalizations of the expressions (36), (37), and (41) from the paper by 

Chandrasekhar (1964). The pulsation equation (43) of his work can be rewritten in our case 
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as 

(57) 

Eqs. (j5ip and (j57p fully describe radial pulsations of superfluid neutron stars. What remains 
to be done is to find the unknown functions 5P and entering these equations. 



4.2 The functions 6P and dji^ 

With the quasineutrality condition, which is valid in a pulsating neutron star, any thermo- 
dynamic function (for a stellar core composed of neutrons, protons, and electrons) can be 
represented as a function of three thermodynamic variables, say, nb, ng, and S (the quadrat- 
ically small dependence of the thermodynamic parameters on w'^^-^W{^k)a is neglected). Since 
the pulsations are assumed to be small, the pressure P{n\,, rie, S) = Pq + 6P and the neu- 
tron chemical potential /in('^b, ne, S) = fino + ^fJ^n can be expanded in the vicinity of their 
equilibrium values, 

Cnbo CrZeo Cdq 

dHn{nw,neO,So) ^ , 9/in(nb0, neO, 5'o) 9;Un(nbO, ^eO, S'o) , . 

= + ^9^^ ^"""^ 7)^ ^^^^ 

Let us find 5nb, 5ne, and SS from the continuity equations for baryons (jH^ . electrons (jTj), 
and entropy (|T3|l . respectively. Writing explicitly the covariant derivative in the metric of 
Eq. (|18p and keeping only terms linear in the perturbations, one can rewrite the continuity 
equation for baryons (jHUj) as 

at or ^ '2 at 2 dr 

1 ^ 1 d 

"^-2 ^'^'^'^W) + 2 d^^ + ""'^ = °- 

The time integration of this equation using Eqs. (^Hj) . and the equality with 

= will yield 

p!^o/2 a pi/o/2 Q 

The expressions for 871^ and 55" can be derived in a similar way: 

5ne = -!^|-(r^neoee--/^), (62) 
5S = -!^|-(.^^oee--/^). (63) 
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Eqs. ()6H ) - ()63|) are generalizations of Chandrasekhar's (1964) equation (50). Note that, Eq. 
()fil|) can be rewritten in a more compact form. By multiplying its left- and right-hand sides 
by fino and using Eqs. dSHI), (I2Z|), dSH), dSS}, and (jnEI), we find 

/^no (^^b = S£ — Tq 6S. (64) 

This is just the second law of thermodynamics © with the quasineutrality (^co = ripo) and 
beta-equilibrium {6fio = 0) conditions valid for an equilibrium star taken into account. 

The substitution of Eqs. into and (jSIJjl gives, after standard transforma- 

tions, 

dP 

SP = -^e-Tii'o^-APo^, (65) 

6Hn = - - /52 yWnO^, (66) 

dr 

with 

$ = -^^{r'^e-'^'^/'), vl> = ^^(rV.oK.^.e-^«/^), (67) 
or n^^or or 

n^o dP{nw,Xeo,Xso) nbo <9/in(?^bo, a;eo, a^so) 
7i = S a ' ^2 = J, , (68) 

^ _ ^ dP{nw,neO, Sp) ^ _ !^ ^/in(^bO, ^eO, 'S'o) 

^ Pq dn^o ' ^ /ino ^'^bO 

where Xeo = '"-eo/'^'bo and Xso = Sq/ui^q. It should be noted that the partial derivatives in Eq. 
(|68|l are taken at constant values of Xeo and Xgo- The new parameter 71 is just an adiabatic 
index of matter that describes pulsations of normal (non-superfluid) stars. When calculating 
the partial derivatives of thermodynamic parameters, one can neglect the temperature effects 
and put So = and Xgo = everywhere. 

Thus, we have found the functions 6fi^ and 6P under the assumption of frozen nuclear 
composition. In this work, all actual calculations of the eigenfrequency spectrum are based 
on this assumption. Still, we would like to make a comment on how one could find these 
functions in the opposite case when uj <^ l/r (when the core is in beta-equilibrium during 
pulsations). The pressure P and the neutron chemical potential fi^ are then functions of 
nb and S only, whereas the electron number density ne(?T-b, S) is derived from the beta- 
equilibrium condition. Using Eqs. (jUT|) and (jU^ . one can write 

6P = ^^|i^5nb + ^%i^55, (70) 
dn^o dSo 

dfln{niyO, So) (9/in(^bO; 'S'o) . „ , . 

OfJ-n = t: OnbH — 6S. (71) 

dribo oSo 

If we now neglect the entropy dependence of thermodynamic parameters (as is justified for 
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the frozen nuclear composition), we will arrive at a qualitatively wrong result, where one of 
the branches of the pulsation spectrum is missing. Indeed, the pulsation equation ()57j) will 
then depend only on the eigenfunction Tq (see Eqs. ESUEH andHOl)- Therefore, the pulsation 
eigenfrequencies can be found just from Eq. alone, independently of Eq. (jHH). [It will 
be shown in the next section that the boundary conditions for the pulsation equations |^ 
andEZIcan also be formulated only in terms of the eigenfunction Tq] The obtained branch 
of the pulsation spectrum practically coincides with the spectrum of a non-superfluid star, 
while the specifically "superfluid" pulsation modes will be lost. To the best of our knowledge, 
such "temperature" pulsation modes have not been discussed previously in the neutron-star 
literature. 



4.3 Boundary conditions 

Pulsation equations (jHH) and dSZj) together with Eqs. (IHn|) . (p^ . (IHD. (IHT)|1 . and en- 

able one to determine the unknown functions ^n, z^, ^ and the frequency spectrum, provided 
that the boundary conditions are known. 

To formulate the boundary conditions, we should specify the model problem to be solved. 
We assume neutrons to be superfluid inside a sphere of circumferential radius Rq with 
Rq ^ -Rcc; where Rcc is the radial coordinate of the crust-core interface. Outside the sphere, 
neutrons are assumed to be normal. The parameters related to the outer (r > Rq) region 
of the star will be marked with the letter "c". On the stellar surface, we have a standard 
boundary condition: 

P,{R + UR)) = ^. (72) 



which can be rewritten as 



5Pc + ^ec 

dr 



0. 



(73) 



J r=R 



Here R is the circumferential radius of an unperturbed star; S,c is the Lagrangian displacement 
of matter in the outer region. In Eq. ()73p we defined Pc{R) = Po{R) + SPc- All derivatives 
with respect to r at the stellar center must be finite, which means that the following limits 
are finite 



lim^/r < oo, 

r— >0 



lim2;,7r < oo. 



(74) 



The other boundary conditions should be formulated at the superfluid-normal interface. 
First, the electron current at the interface must be continuous. It follows then from Eq. ()62|) 
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that the Lagrangian displacement of "normal" particles is continuous, which leads to 

aRo)=URo)- (75) 
In addition, the energy and momentum currents through the interface must also be contin- 
uous. These conditions lead to the following equalities (see Eqs. HOI - 1121 together with the 
expressions IIHI and I^Ujl : 

P(i?o + e(^o)) = Pc{Ro + UR<))), (76) 

[(Po + £o) ^ + /^nO^ioKi 2;4=Ho = [(^0 + £o) ^4=Ro ■ (77) 

With Eq. (f7H|) . the equalities (fTH)) and (|77jl can be written: 

[6P-6P,U^^ = 0, (78) 

/^iO^m Zi \r=Ro = 0. (79) 

Eqs. (f7H|). (fTij) . (f75|). (f7H|) . and ([7^1) cover all boundary conditions that are to be imposed 
on Eqs. ()51|) and (jFFjl in order to find the frequency spectrum for the present neutron star 
model. 



5 SOUND WAVES IN SUPERFLUID MIXTURES 

Before discussing the numerical solutions to the pulsation equations and (|Ff|). let us 
analyze sound waves in superfluid neutron stars. One would expect the numerical solutions 
to resemble a "plane" sound wave when the number of nodes of the eigenfunctions C, and 
Zi is large, so that the wave number is large, k ~ N/R ^ 1/R. Taking into account the 
estimate uj/k ~ u, where u is the sound velocity, we see that the eigenfrequencies of such 
"sound-like" modes must obey the inequality 

u^u/R. (80) 

We now simplify Eqs. (j^T|) and ()57p to the case of short wavelength oscillations. Since the 
characteristic scale R of variation of the equilibrium parameters (marked with the subscript 
"0" ) is much larger than the characteristic scale of the variation of the eigenfunctions, we 
can neglect the spatial derivatives of the "equilibrium" quantities and rewrite the pulsation 
equations as 

/znoe^°--.;^(^n + = (81) 

_eAo-.o^2^^1 ^ (82) 
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Figure 1. Sound velocities iti.2 (in units of c) as a function of temperature T for two models of the nucleon-nucleon potential: BJ 
v6 (solid lines) and Reid v6 (dashed lines). The u\{T) and 1*2 (T) curves are marked as "normal" and "superfluid" , respectively. 
The neutron and proton critical temperatures are indicated by vertical dot-and-dashed lines. The baryon number density is 
ribo = 3no=0.48 fm~^. 



Under the assumption of a frozen nuclear composition, the functions 6P and 5/in are defined 
by Eqs. and as before, but now we have 



dr nbo dr 

The functions ^ and Zi can be presented in the form 



e = eo(r)e 



i{kr—Lut) 



Zi = Zio{r}e 



i{kr—ujt) 



(83) 



(84) 



The derivatives of the slowly varying functions C,o{r) and Zio{r) can be ignored. By sub- 
stituting the expressions ()84|) into Eqs. (jHTj) and (jH^ . one can find from the compatibility 
condition of the resulting set of equations, a biquadratic equation for the local sound velocity 

u = e^^o-'^o^/^^/fc : 



4 I 
yu + 







T^nO^bO 



iPi - 7i - IW) +72-/^2 



H — - P1I2) 

A^nO^bO 



0. 



(85) 
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This equation has two nontrivial solutions for two possible sound velocities (see Andersson 
and Comer 2001a for a similar discussion). The dimensionless parameter y is defined as 

y = n (Y Y -Y Y \ ~ ^^^^ 

f-nO nn -' pp np -' pn J 

At T ^ Ten, we have Ynn, ^np, ^pn 0, hence, y ~ 'n'bo/(/^no ^nn) oo. In this case the 
roots of Eq. (jHK|l are approximately given by 

u, ^ u, ^ {(3211- f3^). (87) 

The first root describes the velocity of perturbations similar to the familiar sound propagat- 
ing through a medium with non-superfluid neutrons. The second root indicates the existence 
of an additional pulsation mode specific to superfluid matter. For the second mode to be 
stable the condition [32ii ^ /?i72 must be fulfilled. The second pulsation mode vanishes at 
T > Ten (^nn = ^np = ^pn = 0), while the velocity of the first mode is still defined by Eq. 
(jS7j). In that case, the first mode is just the usual sound. 

The results of a numerical solution of Eq. ()85j) for matter with baryon number density 
ribo = 3?T,o are presented in Fig. 1 [uq = 0.16 fm~^ is the baryon number density in atomic 
nuclei). In determining these data we used critical temperatures for neutrons and protons 
equal to Ten = 6 x 10^ K and T^p = 5 x 10^ K, respectively, and employed the equation of state 
of Heiselberg and Hjorth- Jensen (1999) to calculate the thermodynamic parameters and their 
derivatives. The velocities ui 2 (in units of c) are plotted as a function of temperature T for 
two models of nucleon-nucleon potential: BJ v6 (solid lines) and Reid v6 (dashed lines). Note 
that, the choice of the model potential determines the entrainment matrix and, hence, 
the matrix Yik (see the paper by Gusakov and Haensel 2005; the microphysics is described 
by Jackson et al. 1982). The U2(T) curves are marked "superfluid" and the Mi(T) curves are 
marked "normal". One can see that the sound velocity mi(T) is practically insensitive to the 
model potential chosen: the solid and dashed lines in the figure coincide. 

The analysis of Fig. 1 shows that the results of a numerical solution of Eq. (jHH|) are 
generally consistent with the above conclusions. We would like to stress that the sound 
velocity Ui does not significantly differ from that calculated from Eq. (jHTj) even at T <^ Ten. 
It is also important that the velocity of the second mode U2 becomes comparable to the 
velocity ui at low temperatures, in contrast to the case of pure helium II. At T < O.STcn the 
velocity U2 rapidly approaches its asymptotic value U2{T = 0). 

Let us discuss briefly sound in beta-equilibrated matter. It is easy to verify that all the 
formulas derived in this section remain valid, provided that the thermodynamic parameters 
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Figure 2. The velocity of tlie second sound U2 (in units of c) in beta-equilibrated matter as a function of temperature T for the 
same models of the nucleon-nucleon potential, neutron and proton superfluidity, and baryon number density tl^q as in Fig. 1. 



and their derivatives are considered as functions of only the baryon number density ribo and 
entropy 5*0- (We recall that the electron number density Ueo is not an independent variable 
in this case; it is determined by the beta-equilibrium condition.) In particular, the adiabatic 
index is now: 71 = (nbo/-Po) 5P(rabo, a;so)/c^^bo- 

Using ribo and 5*0 as independent variables instead of nbo and Xso in the functions 71 and 
72 and expressing the derivatives 9/in(nbo, 5'o)/9?2bo and 9/in(nbo, 5'o)/(9S'o in Eq. (|H3j) with 
the help of the Gibbs-Duhem relation, dP = Sq dT + ribo d/in, we get 

yu'-^— (s^^ + yn^o^ + Soy—\u' + ^^ ( ^ K . ^ \ =0 

(88) 



ibO 



An approximate solution to this equation can be easily found if we keep in mind that we 
always have Ui ^ U2- 



h- dP / Si dT 

V /inO C'^bO V AinO'^bO V oSq 

Again, the first root describes the velocity of sound in non-superfluid matter (the first 
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sound) and the second root describes the velocity of the so-called second sound. It should 
be noted that the first sound cannot, in fact, propagate with the velocity ui defined by Eq. 
(I89|l . because this velocity is so high that no beta-equilibrium can exist in such a wave. If we 
use Eq. flH^ to describe the sound in a one-component hquid, the expression for U2 in the 
non-relativistic limit will coincide with that for the second sound in liquid helium II (see, 
e.g., Khalatnikov 1989). 

The function U2{T) for matter with baryon number density ribo = 3no is shown in Fig. 2. 
We used the same models of superfluidity and the nucleon-nucleon potentials and the same 
equation of state as in the discussion of sound in matter with frozen nuclear composition. 
The speed of sound was calculated numerically. While doing the calculations we used the 
formula C = T dS/ dT, where C is the heat capacity of superfluid matter (an expression for 
C can be found, e.g., in Yakovlev et al. 1999). 

It follows from Eq. (jSHl) and Fig. 2 that the velocity U2 goes to zero at both T = Ten 
and T = 0. However, beta-processes are so suppressed at low temperatures that the second 
sound will not be able to propagate, because matter cannot approach beta-equilibrium on a 
timescale comparable with the pulsation period. Therefore, the second sound can only exist 
in a range of temperatures near T < Ten- 

To conclude, three types of sound waves can exist in superfluid npe-matter. The speed 
of two of them is so high that they propagate in matter with a frozen nuclear composition, 
while the waves of the third type can exist only in beta-equilibrated matter at temperatures 
in the vicinity of the neutron critical temperature T^n. 



6 RESULTS FOR RADIAL PULSATIONS 

Let us now discuss the solutions to the pulsation equations and ^7\i . We have integrated 
the equations in a standard way, using the Runge-Kutta method. We employed the equation 
of state of Negele and Vautherin (1973) in the stellar crust and that of Heiselberg and Hjorth- 
Jensen (1999) in the core. The latter is a convenient analytical approximation to the equation 
of state proposed by Akmal and Pandharipande (1997). For this equation of state, the most 
massive stable neutron star has central density pc = 2.76 x 10^^ g cm~^, circumferential 
radius R = 10.3 km, and gravitational mass M = Mmax = 1.92 M©. The powerful direct 
Urea process of neutrino emission is open in the core of a star of mass M > 1.83 Mq. When 
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calculating the matrix Yik we have used the model BJ v6 of nucleon-nucleon potential (see 
Section 5). 

For illustration, we consider a neutron star model with mass M = IAMq {R = 12.17 
km, pc = 9.26 X 10^^ g cm~^). For such a star the crust-core interface is at -Rcc = 10.88 km. 
The frequencies of the first three modes of radial pulsations of a non-superfluid star with 
this mass are uoi = 1.703 x 10^ s-\ 002 = 4.081 x 10^ s-\ and uj^ = 5.732 x 10^ s-\ 

To reduce the number of factors affecting the pulsation spectrum, we consider a simplified 
superfluidity model in which the critical red-shifted temperatures of nucleons do not vary 
with the density and are equal to = T^ne""^'^ = 6 x 10^ K and T~ = T^pe^^o/^ ^ 
5 X 10^ K. Consequently, superfiuid matter is contained in the stellar core: Rq = Rcc- This 
means that the boundary at Rq is "attached" to matter and, for example, is independent of 
temperature variations. (Note that, in the more general case of density dependent profiles 
of critical temperatures, the superfiuid-normal boundary can depend on T and temperature 
perturbations.) Numerical tests have shown that the approximation of critical temperatures 
T^p as constant throughout the core describes reality well if these temperatures smoothly 
depend on the density. This is consistent with the predictions of some microscopic models 
of nucleon pairing known in literature (see, e.g., Yakovlev et al. 1999). 

Fig. 3 shows the dependence of the pulsation eigenfrequencies u on the red-shifted tem- 
perature T°° = Toe^"/^ (we recall that the superfiuid core is isothermal, in accordance with 
Eq. EJ. The vertical dot-and-dashed line indicates the neutron critical temperature T^. 
The horizontal dotted lines show the first three eigenfrequencies uji, uj2, and for a non- 
superfiuid star. No attempt to determine the spectrum in the shaded region was made. At 
T°° > T^^, the star pulsates as a normal fluid (no matter whether the protons are paired or 
not). Hence, the spectrum contains only normal, temperature-independent pulsation modes 
(the first three modes I, II, and III are shown as solid lines). At T°° < O.IT^^, a pulsating star 
can be described in the zero-temperature approximation. The spectrum of a cold superfiuid 
star is doubled, as compared with that of a normal star (see Comer et al. 1999). In addition 
to "normal" pulsation modes, whose eigenfrequencies are close to those for a non-superfiuid 
star (solid lines), the spectrum contains specific "superfiuid" modes (dashed lines). Note 
that, the first "superfiuid" mode is quite different from the "normal" one but the second 
and third "superfiuid" modes are already sufficiently close to their "normal" counterparts 
(see Fig. 3). 

As the temperature increases, starting from approximately T°° ~ 10* K, the frequency 
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Figure 3. The pulsation eigenfrequencies uj (in units of ljq = c/i?) of a neutron star as a function of tlie red-shifted core 
temperature T°°. The neutron critical temperature is indicated by the vertical dot-and-dashed line; the horizontal dotted 
lines changing into solid lines at T°° > indicate the first three eigenfrequencies (I, II, III) of a non-superfiuid star. No 
spectrum was plotted in the shaded region. The dashed curves correspond to "superfluid" modes at T°° <^ and the solid 
curves correspond to "normal" modes at T°° <C (see text). 

of each mode begins to decrease. When a mode reaches one of the horizontal dotted hnes, 
it changes behavior and becomes temperature-independent, imitating the behavior of one 
of the non-superfluid modes. As the temperature rises further, the frequency of the higher 
mode approaches that of the mode in question, which in turn begins to decrease again (see 
avoiding crossings in Fig. 3). As a result, the two different modes of the spectrum will never 
intersect. One can conclude that a given mode may behave either as "superfluid" or "normal" 
with increasing temperature. 

The behavior of the frequency spectrum at temperatures close to is of particular 
interest. It is clear from Fig. 3 that the frequency of any mode goes to zero at T°° = T^. 
This is not surprising if we keep in mind that high order pulsation modes represent sound-like 
waves (see Section 5), and that the frequency of the "superfluid" sound also goes to zero at 
the transition point into the superfluid state (Fig. 1). It might seem that the spectrum does 
not contain eigenfrequencies of non-superfiuid stars at the transition point when all neutrons 
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in the star are normal. However, this is not the case. The point is that at — ^ T^, the 
number of modes with frequencies in any given interval, say, [0, ui], becomes infinitely large. 
As a result, at any temperature T°° and any eigenfrequency of a normal star, there is a 
mode which is temporarily "normal-like", i.e., it has the same frequency as in the normal 
fluid. 

Since the temperature of neutron stars changes with time, it would be interesting to 
discuss how the pulsation frequencies vary with time. Suppose that the pulsation energy 
is much lower than the thermal energy. We can then neglect the star heating due to the 
conversion of pulsation energy into heat (see Gusakov et al. 2005 for details). The star will 
cool down, and to determine the dependence of the internal temperature T°° on time t one 
should use the cooling theory of superfluid neutron stars (see, e.g., Yakovlev et al. 1999, 
Yakovlev and Pethick 2004). 

Since the direct Urea process is forbidden for the chosen neutron star model, the main 
coohng mechanisms (at T°° < T^) will be neutrino emission due to Cooper pairing of neu- 
trons, neutron-neutron bremsstrahlung, and the photon emission from the stellar surface. 
One can easily find the function T°°(t) by solving the thermal balance equation (see, e.g., 
Yakovlev et al. 1999) under the assumption that the stellar core is isothermal. If the depen- 
dencies c<j(T°°) and T°°{t) are known, it is possible to plot the frequency spectrum o; as a 
function of time t (Fig. 4). Here the time (in units of 10^ years) is counted from the moment 
of neutron superfluidity onset (at T°° — T^) . 

The analysis of Fig. 4 shows a significant change in the pulsation spectrum for 20 years 
after superfluidity turns on. This is associated with the highly efficient Cooper pairing neu- 
trino emission (the detailed discussion of this process and its influence on the neutron star 
cooling is given by Gusakov et al. 2004). For example, the frequency of the third "superfluid" 
mode changes during this period of time from to the eigenfrequency CU2 of a non-superfluid 
star. The Cooper pairing neutrino emission process quickly becomes weaker with time, the 
cooling slows down and the variation in u(t) becomes smoother. We would like to emphasize 
that the fast change of the pulsation frequencies for the first few dozens of years is due to the 
high critical temperature of neutrons, = 6 x 10^ K. We could make the Lj{t) dependence 
less dramatic by choosing lower critical temperatures. 
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Figure 4. The dependence of the pulsation spectrum of a superfluid neutron star on time t, counted from the moment of 
neutron superfluidity onset (at T°° = T^). The time is in units of 10^ years. Notations are the same as in Fig. 3. 



7 SUMMARY 

The aim of the present study was to analyze radial pulsations of superfluid neutron stars 
at finite core temperatures. We used the equations for one-component superfluid hydro- 
dynamics suggested by Son (2001), rewritten in terms of more convenient variables, and 
extended to the case of superfluid mixtures in General Relativity. A simple model of npe- 
matter was employed to show that a necessary condition for a star to be at hydrostatic and 
beta-equilibrium is constancy of the red-shifted temperature in the region of the star where 
the neutrons are superfluid: Te'^"/^ = constant. Proton superfluidity does not impose any 
restrictions on the temperature, because protons are "coupled" with normal electrons by 
electromagnetic forces and behave as a normal fluid, no matter whether they are superfluid 
or not. 

The hydrodynamics of superfluid mixtures was applied to investigate radial pulsations 
of neutron stars. It was assumed that the crust is non-superfluid, and neutrons and protons 
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have red-shifted critical temperatures, which are constant throughout the core. The set of 
equations we have derived describes radial pulsations of superfluid stellar matter. 

We have found the short wavelength solutions to this set of equations, representing sound 
waves in superfluid neutron star matter. The dependence of the speed of sound on the 
stellar temperature was examined in two limiting pulsation regimes: (1) in beta-equilibrated 
pulsating matter and (2) in pulsating matter with frozen nuclear composition. It was shown 
that three different kinds of sound waves may in principle exist, two of them propagate in the 
matter with frozen nuclear composition and one can exist only in beta-equilibrium. While 
the speeds of the former sound waves are comparable to each other (see Fig. 1) and to the 
speed of sound in the usual non-superfluid matter, the speed of the latter is 4-5 orders of 
magnitude lower (see Fig. 2); it can be excited only at temperatures T close to Ten- 

Generally, the pulsation equations were solved numerically, and the results show that 
the finite internal temperatures strongly affect the pulsation spectrum in the range of T ~ 
(0.1 — l)Tcn (see Fig. 3). The frequency of any pulsation mode in this range decreases with 
increasing temperature. However, when the mode reaches one of the eigenfrequencies of 
a non-pulsating star, it becomes temperature independent for a while. One may say that 
it begins to mimic the behavior of a non-superfluid mode. At T — > Ten, all superfluid 
eigenfrequencies tend to zero. At T < O.lTcn, the pulsation spectrum is similar to that 
calculated in the zero temperature approximation. 

In addition to the analysis of the temperature dependence of the pulsation spectrum, we 
discuss the temporal evolution of the eigenfrequencies during the star cooling (Fig. 4). In our 
analysis, we use the standard cooling theory of superfluid neutron stars (see, e.g., Yakovlev 
et al. 1999). The calculation shows that essential changes (within the present model) in 
the pulsation eigenfrequencies occur for the first 20 years following the moment of neutron 
superfiuidity onset. This rather short (for the cooling theory) period of time is associated 
with the fast cooling due to the effective Cooper pairing neutrino emission process. It will 
be even shorter if the powerful direct Urea process operates in the stellar core. 

The consideration of the problem presented here is based on a simplified model. In partic- 
ular, we discuss only the simplest case of radial pulsations and assume critical temperatures 
of nucleons that are constant throughout the core. However, it would be important (and 
interesting) to understand how finite internal temperatures affect the frequency spectrum of 
non-radial pulsations and how the results would change if we analyzed more realistic density 
profiles for the critical temperatures. Finally, in a more realistic approach one should take 
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into account ^5*0 neutron pairing in the stellar crust and more accurately treat the physics of 
the crust, especially if one deals with pulsation modes localized in the outer layers of the star. 
In spite of the considerable simplification of the problem discussed in this paper, we conclude 
that finite internal temperatures significantly affect the pulsation spectrum of not too cold 
superfluid neutron stars. Moreover, the pulsation frequencies can change dramatically for a 
period of several dozens of years, an effect that may potentially be observable. 
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